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TECINICAL NOTE 2891 

FACTORS AFFECTING LAMINAR BODlIDARr lAIER 
MEASDRMENTS US A SUPERSONIC STREAM ' 
By Rotert E. Blue and George M. 'Lavr 


SDMMARY 

The ohserved discrepancy at supersonic speeds Between theoretical 
and apparent experimental average flat plate friction-drag coefficients 
calculated from Boundary layer tgtal-pressure surveys was investigated. 
Effects of the total-pressure proBe, heat transfer throu^ the leading- 
edge region^ change in leading-edge radius and strength of the leading- 
edge \jave^ possiBle early traneitlcn to turBulent flow or "Bursts of 
turBulence, " and the slight streamwise pressure gradient Inherent in 
flat plate flow were investigated for plates with very sharp leading 
edges. 

Only one of these factors, the effect of the total-pressure proBe, 
was found to Be significant. Total-pressure proBes of different tip 
heights, when placed in laminar Boundary layers developing under 
identical conditions, were found to yield different values of friction- 
drag coefficient. Extrapolation, of these measurements Indicates that 
a proBe of vanishing size would yield the theoretically predicted 
values of average flat plate friction-drag coefficients. A correlation 
descriBlng the relation Between the friction-drag discrepancy and the 
probe tip height is presented. 


INTRODUCTION 

Most available experimental values of average friction-drag 
coefficient for laminar flat plate flow have Been obtained By total- 
pressure surveys of the Boundary layer. Results reported in refer- 
ences 1 to 3 -consistently indicate diag coefficients that are in the 
range of 10 to 100 percent higher than the values predicted By com- 
pressible .laminar Boundary layer theory. Other measurements made along 
the wall of a supersonic tunnel Behind a Boundary layer .remo-val scoop 
(ref. 4) gave values that were even higher than those found in the 
investigations of references 1 to 3 at ccanparaBle conditions of Reyaolds 
number . 



2 


MCA TN 2891 


The wide range of Trarlahillty of the experimental results tends to 
rule out the possihility of attributing the discrepancy between theory 
and e35)erlment to a deficiency in the theory alone. Furthermore, local 
frlcticai coefficients determined from force measurements on a small 
element set in a plate surface show excellent agreement with boundary 
layer theory (ref. 5). The theory considers boundary layer growth 
on an infinitely thin flat plate starting at zero thickness at the 
leading edge and developing into a flow vrith zero pressure gradient. 
Experimentally, this situation can only be approximated . The precision 
to which any part or all of the physical situation hypothesized in the 
theory is reproduced may Influence the size of the measured discrepancy. 

An analysis of local friction coefficient based on the measured 
boundary layer Mach number distribution led the authors of reference 3 
to conclude that the boundary laysr air suffers a momentum loss near 
the plate leading edge. A consideration of the pressure distribution 
over the leading- edge region, based on the blunt-body theory of 
reference 6, shows that such a loss should occur and that it will be 
directly proportional to the radius of the edge. The discrepancies 
in average friction coefficient reported in reference 3 are at most 

3^ times larger than those shown in reference 1. The leading-edge 

radii, however, differed by a factor of 10. Conseq.uently, factors 
other than the leading-edge pressure drag must be involved in the 
observed discrepEuncy. 

The measured momentum discrepancy could be attributed to any of 
a number of causes, for example, heat transfer throu^ the thin 
leading-edge region of the plate, momentary transition to turbulence, 
that is, "bursts of turbulence" (ref. 7), or an instrument error in the 
form of a probe interference effect or a misinterpretation of the pres- 
sure reading. It is the pirrpose of this report to present the findings 
of an investigation of these and other effects vThich mi^t account for 
the difference between theory and experiment. 

This investigation consisted principally of an eaperimental study 
of the boundary layer on a flat plate and an evaluation of the 
measuring technique. Total-pressure surveys of the boundary layer 
were made at Ifech numbers of 2 and 3 over a range of Eeynolds number 
at several distances from the plate leading edge. Several leading- 
edge radii -tfere used on the flat plates and a systematic variation 
of total-pressure probe tip height was made. These measurements vrere 
made in two different test facilities at the HACA Lewis laboratory. 

The authors are indebted to Mr. John W. Schwartzenberg for TnalHng 
many of the experimental measurements. 
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APPARATDS AH3) PEOCEDUEE 

The ‘boundary layer measurements -were jpade on flat-plate models In 
two supersonic wind . tunnels . Tunnel 1 had a test section of 
3.84 "by 10 inches and a nominal free-stream Lfech number of 2. The 
test section of tunnel 2 measured 4 hy 4 Inches, and its nominal 
free-stream Mach number iras 3. Air to both tunnels tras supplied at 
^ a maximum staguation pressure of 50 inches of mercury absolute. The 

^ supply air temperature (ahead of heaters) was -20° F for tunnel 1 and 

N -70° F for tunnel 2, thus insuring a dew point sufficiently low to 

obtain condensation-free flow. The tunnels eshausted to a pressure of 
about 3 inches of mercury absolute. 

All total pressures In the boundary layer were measured on a 
mercury well-type manometer. Static pressures were measured on butyl 
phthalate manometers irtth hi^ -vacuum reference pressures. Temperature 
measurements were Indicated on a self-balancing potentiometer. 

The teat plate used In tunnel 1 was 6 Inches long, 3.84 inches 
wide, and 0.25 inch thick. The plate used In tunnel 2 measured 
8 by 4 by 0.25 inch. The top (test) surfaces of both plates were flat 
and the surface roughness did not exceed a -value of 8 microlnches rma. 
Eubber seals were pro-vlded be-tween the plates and the tunnel side walls 
to prevent flow be-tween the top and bottom surfaces. 

The plate used in tunnel 1 had a leading-edge radius of curvature 
of less than 0.0005 Inch, as determined on a comparator from a section 
of a molding plaster cast of the leading-edge region (a 50X magnifica- 
tion of the leading-edge region shovred that the nbse of the plate was 
nearly semicircular, so that definition of leading-edge radius tmB 
justified). The bottom surface of the plate was chamfered to an angle 
of 12° in the leading- and traillng-edge regions. 

In order to determine the pressure gradient and the tempera-bure 
along the plate, 15 s-batic-pressure taps and 6 thermocouples were 
located in the plate surface, as sho-wn in figure 1. The plate -was 
supported by pivots and a sting in order to permit changes in angle 
of attack. 

The radius of curvature of the leading edge of the plate used in 
tunnel 2 was less than 0.0005 inch for several tests, but a larger 
radius of cur-vature, of the order of 0.001 to 0.002 Inch, was used in 
most of the tests. The bottom surface was chamfered to an angle of 
10.6°. Static-pressure taps were Installed in the plate surface 
2.38 and 3.38 inches from the plate leading edge and 0.25 inch from 
the plate center line. 
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Details of tlie total-pressure protes used to survey the ‘boundary 
layer are shown In figure 2. A range of prohe tip heights from 
0.003 to 0.016 inch is represented. Distance from the plate surface 
to the bottom of the prohe tip ■was measured -tvlth a micrometer and is 
considered accurate to 0.0005 inch. Contact of the probe with the 
plate surface ims Indicated electrically (see ref. l) . 

Only probe 3 (see table, fig. 2) was used in tunnel 1. Total- 
pressure surveys of the boundary laysr were made with the probe tip 

located 1^, and 2^ inches from the plate leading edge with the plate 

at zero angle of attack, and with the probe tip located 2^ inches from 

the leading edge at an angle of attack of -5°. For each probe location, 
the inlet stagnation pressure was varied from about 15 to 40 inches of 
mercury absolute. The. inlet air was heated to a stagnation temperature 
ranging bebtreen 80° and 90° F. 

Probes 1, 2, 4, and 5 wore used in tunnel 2, and boundary layer 

3 3 

surveys were made with the probe tip located 2-g’ and 3g inches from the 

plate leading edge. The inlet stagnation pressure was varied from 
20 to 40 Inches of mercury absolute, while the inlet temperature ■jvas 
held at 102° F. 

Plate static pressures and surface temperatures were recorded 
before and after each boundaa:^ layer survey. All surveys w^re started 
at the plate surface. 


DATA REDUCTION 

Calculation of momentum lose in the boundary layer req.uiree the 
evaluation of mass flow pu and momentum pu^ of the boundary layer 
air. In order to reduce the calculation steps req,uired to obtain 
these q.uantities from the experimentally determined values of static 
and pitot pressure, the following dimensionless parameters were used 
(ref. 8): 



velocity parameter 

p 

mass-flow parameter 


where T is the total temperature. (All symbols used in this report 
are defined in appendix A.) The mass-flow parameter was introduced in 
reciprocal form as a static-pressure parameter in reference 8. These 
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paxameters can "be slioi'in 110 , 1)0 uniq.uely determined t)y tlie ratio of static 
to pitot pressure p/H. The calculations of and aP = pu^/p from 

p/H were made hy machines in the manner described in appendix B. 
Typical laminar boundary layer profiles of the Telocity, mass-flow, and 
momentum parameters are shown in figures 3 to 5 for seTeral Reynolds 
numbers. These data represent all the boundary layer surreys made with 
probe 4. Duplicate surreys under identical conditions are distinguished 
by open and solid points. There was a sli^t lack of reproducibility, 
as is apparent in these figures, for all the boundary layer surreys m^e. 

The mcjnentum thickness 6 expressed in terms of the paioameters a 
and p is 



With the conrentional assumption that static pressure is constant throu^ 
the boundary layer, the ratio p/p]_ is 1. The distribution of total 

temperature is also required to eraluate the equation for 6. Because 
no temperature surreys were made, the simplifying assunption of a constant 
total temperature was used in all calculations. The, effect of this 
assumption on the calculation of 9 is discussed in the section 
Calculation Errors. 


The total-pressure reeidings, assumed to exist at the geometric 
center of the tube, were taken at equally spaced distances from the plate 
surface. It was therefore possible to integrate the equation for 0 
numerically ’tilth the a emd p calculated from the original data. 

Values for the integrals of p and op were obtained for erery other 
point of the original data using Stirling's rule. 


A "running ralue" of momentum thickness 6j tiras defined by replacing 
the free-stream quantities subscripted 1 in equation (l) with ralues 
measured at the ordinate y, denoted by a subscript y in equation (2), 
and carrying the integration out to the limit y. 


0 


J 



= rr- 


1 


Py py 

P dy - op dy 

0 ^^7 Jo 



The ratios of pressure and temperature were eliminated in equation (2) . 
by the assumptions of constant static pressure and total temperature 
through the boundary layer. 

The method just described of calculating a profile of 0y rather 
than predeteimlnlng the free-stream relocity and temperature or Mach 
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numljer as ■was done in references 1 to 3 shows th.e dependence of 6 on 
the choice of free-stream conditions. The profiles of dj shovm 

in figure 6 attain a peak -mlue before reaching the edge of the boundary 
lajBr and then decrease, approaching a constant in the stream far from 
the plate surface. Yal'u.es of 0 ueed to calculate friction-drag 
coefficients were chosen at the point where the velocity gradient first 
became zero. Average flat plate friction coefficients were calculated 
from the Earm^ momentum eq.uation 


2 


^ + 20 
dx 




(3) 


Integration of over a length of boundary layer run x yields the 

average friction-drag coefficient. For the case of a zero pressure 
gradient the momentum eq.uation (3) may be readily Integrated , yielding 
Cj = 20/x for the average friction-drag coefficient. 


Experimental Errors 

Static pressure . - Pressures measured ■vrith a differential manometer 
board using butyl phthalate as a working fluid can be read to an equiv- 
alent accuracy of better than ±0.005 inch of mercury. For the range of 
static pressures in the test sections of the 'vrlnd tunnels used for these 
measurements, this accuracy corresponded to a determination of the static 
pressure -within ±0.8' percent or less. 

A calculation -was made of the effect of an error in static pressure 
at a Mach number of 3 on the quantity -^Ee using the theory of 
reference 9. "Assumed values of Reynolds number, T-^, p, and x which 

were comparable irtth those of the experimental investigation were used 
to calculate a theoretical Mach number distribution through the boundary 
layer. A profile of p/S was then found from the Eayleigh equation and 
the Mach number distribution. Multiplication of the p/h curve by a 
new value of static pressure divided by the original assumed -value gave 
a new distribution of p/s against y which -was used to calculate a 
new value of 0 in the manner described for a measured profile. A 
Reynolds number corresponding to the erroneous static pressure was also 
calculated. The results of this error calculation are presented in 
figure 7. 

A similar error calculation was made using experimentally determined 
distributions of p/H rather than the theoretical profile of reference 4. 
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20 , 

Curves sliOTfiug tlie error in — ■yJ'Re^ caused "by an error In static 

pressure using measured profiles for two different prole sizes are also 
plotted in figure 7. The shape of the boundary layer profile apparently 
has a large effect oq the eiror. The smallest probe used (probe l), 
which measured a profile closest to that of the theory, shows the same 

trend of an eiror in ifith static-pressure error as does the 

error calculation mhde from the theoretical profile. The largest 
probe (probe 5), however, shows a reverse trend. By either of these 
calculations a reading error of 1 percent in static phessure will cause 

an error in — -y/B^ of 1 percent or less. 

Total-pressure recording . - The maximum^ possible error in the 
pressure measured by the pitot tube would be ±0.02 inch of mercury. 

This error would be randomly scattered through the readings taken 
during a boundary- layer survey. No attempt was made to calculate the 
error in 6 due to such a random scatter. It was noted, hcrtrever, that 
an error in the barometer reading, from which the total-pressure reading 
■v/as subtracted, could be as large as l/lO inch of mercury and still affect. 
6 by less than l/2 percent. In the absence of a method for specifically 
evaluating the effect of a random reading error, this result was assumed 
to Indicate that an error in reading either the barometer or the total 
pressure could be neglected. 

Turbulence level . - As sho™ in figure 1 of reference 10, tunnel 1 
has a set of turbulence screens in the surge tank -where the velocity is 
of the order of 8 feet per second. These screens are followed by a 
smooth contraction of 19 to 1 to the tunnel throat. The inlet to 
tunnel 2 consists of a 1-foot-diameter pipe connecting a 2-foot-diameter 
surge tank to a bellmouth. A contraction of 27 to 1 to the tunnel throat 
follows . A dirt filter made from 16 thicknesses of cheesecloth was 
placed across the connecting pipe at a point -vdiere the air velocity was 
about 70 feet per second. 

A conclusion cannot be made regarding the levels of turbulence in 
the two tunnels from the measurements taken. The extreme difference 
in design, however, can be assumed to have caused substantially 
different levels . The agreement of the results obtained in these 
tunnels thus implies that turbulence level has no noticeable effect 
on 0, at least -vdiere the level of turbulence is reasonably low. 

Total -pressure probe . - The possibility exists that a small total- 
pressure probe of the type used for boundary layer surveys fconsistently 
reads a value of total pressure that is in error, even out in the free 
stream. In order to determine whether such an error exists, a 
conventional blunt nose total -pressure probe, l/8 inch in diameter. 


was 
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positioned irith its center line l/2 incli above the plate surface. The 
readings of this prohe in comparison \rith the readings of the various 
houndary layer prohes "when they were located at the same position in 
the flow are shown in figure 8 for the five inlet pressures used in the 
number 2 tunnel measurements . The solid lines in the figure represent 
the readings of the l/8 -inch-diameter prohe. 

The boundary layer prohe readings tend to scatter about a mean 
value slightly Icrtrer than that read by the large probe. The largest 
difference in total pressure measured by the two types of probe was 
2 percent. This correspcnde to approximately an 0.8 percent error in 
Mach number, a, p, or 0. 

Yaw of total-pressure probe . - Each time a total-pressure probe 
is Installed in the trlnd tunnel, it is necessary to realine the probe 
ifith the stream direction. This allnement was made both visually and 
by adjusting for the TnaximiiTn total-pressure reading. Since this type 
of boundary layer total-pressure probe is relatively insensitive to 
angles of yaw in the free stream, it was difficult to be certain thQ 
probe was pointed precisely into the stream. Because the possibility 
that an angle of yaw mi^t cause a more serious error when the probe 
is Immersed in the velocity gradient of the boundary layer tos 
recognized, a total-pressure survey was made for which the probe Tras 
deliberately mlsalined 3°. A running plot of momentum thictoiess 
against y for survey at 3° yaw is compared with the profile tahen 
with the probe alined with the stream in figure 9. The two curves 
agree within less than the experimental reproducibility found for most 
of the boundary- layer surveys (see solid and open points of fig. 6) . 

The actual error in alining the total-pressure probe was less than 
the 3° used in this 'error measurement for the effect of yaw. 


Calculation Errors 

Total -temperature assumption . - As previously stated it is necessary 
to assume some fom of a total-temperature distribution through the 
boundary layer in order to calculate 0, All the calculations for this 
report assumed a constant total temperature. In order to evaluate the 
effect of this assumption on 0 the theoretical Mach number distributions 
in the boundary layer as calculated from reference 9 were used to obtain 
the values of a and 8 necessary to Integrate eq.uation (2) . The 
difference between calculated 0 for a constant total temperature and 
theoretical 0 is plotted in figure 10 In terms of the percentage error 
In 0 against Mach number. A similar calculation was made iflth the 
assumption of a linear variation of total temperature between the 
theoretically predicted wall temperature and the free-stream total tem- 
perature. Errors due to this foim of an assumed temperature distribution 
are also shown in figure 10. It is apparent in this case that the 
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simplest temperature assun^jtion gives tlie test results. At a Maeli 
numter of 2 tlie error is nearly zero and it is still un4er 1 percent 
at a Ifech numter of 3. 

Calculationa near surface . - The typical "boundary layer profiles 
of figures 3 and showing the q.uantities ^ and plotted 

against j, shovr that the first fmr points measured near the plate 
surface do not fall on a smooth, curve approaching the origin. (This 
has "been tentatively attributed to Interference "betrreen the probe and 
the plate surface.) Under the assumption that the velocity must go 
to zero at ,the plate surface, fairings neglecting these few apparently 
erroneous points were made. IVom these it was determined that either 
the fourth or the sixth points would lie on a straight line through 
the origin that fairs into the 3 curve ^d on a parabola through 
the origin that fairs into the <xp curve. As explained in appendix B, 
this stral^t line and parabolic variation of p and op was used in 
the integration of these q.uantities out to the fourth or sixth points. 
Beyond these points, Stirling's rule tra.s used to Integrate between the 
measured data points. This method of fairing is compatible "vrlth the 
theoretically predicted change In velocity and density near the plate 
surface. Since the velocity near the surface changes linearly while the 
density is nearly constant for the condition of zero heat transfer, the 
product pu will vary almost linearly. The pu^ term, being the product 
of ti7o linearly varying q.uantities, behaves as a parabola. Several 
calculations of 6 steirtlng at both the fourth and sixth points of a 
profile were made to see what difference the choice of starting point 
would have on 0. It was found to affect the calculated value of 6 
by about 0.5 percent. 

Probe zero location . - An error in setting the zero location of 
the probe Kirill cause an error in the Integrals of p and op as a 
result of the method of computation used for the first few points. This 
can be seen from the following sketch of p against y. The points 
which are eliminated from the. calculation of 9 by the fairing method 
are numbered 1, 2, and 3. 
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If the zero setting is off hy an amount Ay, the additional cross-hatched 
area hounded hy the straight line to the origin and hy the broken line, 
which is the original curye displaced hy the error Ay, will he added to 
the Integral of p. It is eq,ual to 84 Similarly, for op the 


integral tjIII he in error hy 


^ 4 :^ 4 = 

o 


the error in 0 will he 


Ay. The theoretical yariatlon of P and op computed from 


reference 9 for an awerage Eeynolds numher and Ifech number of these 
experiments was used to calculate the percentage error in 0 due to an 
error in y at zero for fairings starting at seyeral distances from the 
wall. The results eire plotted in figure 11. The error in locating the 
probe zero position for the present measurements was less than 
iP.0005 inch, which amounts to an error of ±1 or 2 percent in 0 for the 
range of starting points In y that ms used. 


The following table summarizes the maximum -values of all possible 
errors in 0 : 
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Cause of error 

Eiror 

(percent) 

Static pressure 

±1.0 

Total-pressure probe 

-0.8 

Total-temperature assumption 

-0.7 

Probe zero location 

±1.6 

Total error In 9 

+1.1 

-4.1 


2g 

It "was necessary to take the error in -y/Re^ caused hy an error In 

static pressure to he the same as the error In 6 for this estimate. 
The reproducihility seen in figure 6 is -within the range of this max- 
imum to-bal error. 


EEBUIiCS AJND DISCUSSIOU 

Average flat-plate friction-drag coefficients 2e/x -were computed 
for all measured total-pressure profiles and are plotted against local 
stream Eeynolds number Ee^ on logarithmic coordinates in figure 12. 

A theoretical -variation of Ze/x against Ee^- obtained from reference 9 
is also shown in figure 12 and applies for both = 2 and M-l = 3 

under the inlet temperature conditions of these tests. The measured 
values of ZO/x are from 5 to 65 percent higher than those predicted by 
theory. These high -values of friction-drag coefficient might be 
attributed to transition to turbulent flow, to the sli^t longitudinal 
pressure gradient inherent to flow over a flat plate, to a momentum loss 
originating at the leading edge of the plate, to heat transfer -throu^ 
the leading-edge region, or to an effect of -fche total-pressure probe. 

Transition from l^lnar to turbulent flow . - In order to determine 
the exbent of the laminar boundary layer, the inlet stagnation pressure 
-was -varied while the to-bal-pressure probe was held in a fixed position 
in the boundary layer. As long as the boundary layer isl^lnar, points 
on a plot of op against the similarity parameter jtilu.^/v^x should fall 

near a universal curve regardless of the value of y. Any large de-vlation 
from this curve, which occurs at different values of op depending on 
the distance y from the plate surface, indicates transition from laminar 
to turbulent flovr. Typical transition curves ob-balned in tunnel 2 at 

3 

3: = 3g inches are sho^m in figure 13 together -with a conventional op 

profile. Transition occurs at Ee2-~10®. All results presented in this 

report were obtained at Eeynolds numbers less than the transition 
Eeynolds number as indicated by curves similar to figure 13. 
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The possibility that "bursts of turbulence" (ref. 7) are respon- 
sible for the momentum discrepancy while the mean value of the pro- 
files did not exhibit transitional tendencies was also investigated. 

A low-range pressure transducer X'jas connected to the total -pressure 
probe with the shortest possible length of tubing. The output of the 
transducer ■tra.s fed throu^ a preamplifier to an oscilloscope. The 
system could respond to pressure fluctuatione as fast as 130 cycles 
per second. Ho pressure fluctuations, however, were detected at any 
point in that portion of the boundary layer where profiles were 
obtained. It therefore appears unlikely that momentary transition to 
turbulent flon^ existed in the Eeynolds number range of these experiments. 

Effect of pressure gradient . - Pressure distributions obtained on 
the flat plate In tunnel 1 are presented in figure 14. The Initial 
expansion and f oUovring compression near the leadlog edge are character- 
istic of the flat plate flow and are not /due to a tunnel disturbance. 

This fact was established by moving the plate upstream and downstream 
from the teat position. The con^jresslon near the trailing edge, 
however, is caused by a tunnel disturbance, but it occurs 2 Inches 
downstream of the last test station. 

Althou^ the actual pressure changes are minor, they occur over a 
short distance x, so that the pressure gradients are q.uite large. 
According to theory, a pressure gradient affects not only the momentum 
thickness but also the drag coefficient, which in the presence of a 
pressure gradient is not eq,ual to 2e/x, as can be seen from the Karm^ 
momentum eq.tiation (eq.. (3)). 

Values of momentum thickness and drag coefficient for the pressure 
distributions presented in figure 14 were calculated by the method of 
reference 11 emd compared with a zero pressure gradient solution. 

Results of this comparison sho'vred that the effect of the pressure 
gradient on the momentum thickness was sli^t (1.6 percent), and that 
the error Introduced by letting the drag coefficient eq.ual 2©/x was 
also very small (2,4 percent). 

Effect of leading edge . - If the discrepancy between theoretical 
friction-drag coefficients and measured drag coefficients is assumed 
to be due to an initial momentum defect which has been falsely attri- 
buted to friction, then a change in leading-edge geometry or a change 
in the strength of the leading-edge wave would be expected to change 
the magnitude of this discrepancy. The following tests were made to 
check this assumption: 

The effect of the strength of the leading-edge wave was investigated 
In tunnel 1 by changing the eingle of attack of the plate and thereby 
changing the strength of the shock trave. Because of tunnel choking 
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characteristics, tlie plate angle could he -varied only he-bween 0° and -5°. 

(A negati-ve angle of at-fcack yields a compression at the leading edge . ) 

* / 

A comparison of "velocity profiles at angles of at-back of 0° and -5° 
shows that the profiles were unaffected hy this change in angle of at-back 
(f ig. 15) . The ahscissa in figure 15 is the similarity parameter 
y jJu-iAiZ. Il’ee-stream conditions for this comparison -were selected behind 

the leading-edge shock -waye. Drag coefficients, which are plotted in 
figure 12 also show no change "with angle of attack. 

It could he argued, however, that a detached -wave stands ahead of 
the leading edge, and that the momentum defect caused hy the pressure 
drag of the leading edge is responsible for the high measured drag 
coefficients. On this basis it is not sturprising that the drag 
coefficients show no appreciable change -with angle of at-back, and it 
would be espected that the leading-edge geometry is of prime importance. 

The plate in tunnel 2 \jslb therefore tested -with .an exceedingly 
sharp leahlng edge (r < 0.0005) and -wlti. a somewhat duller leading edge 
(0.001 < r < 0.002). Telocity profiles obtained "vrith both leading-edge 
configurations are presented in figure 16, and corresponding plots of 
against y are shown in figure 17. It is apparent that the 

difference be-tween profiles obtained with the dull and shaip leading edges 
is of the same order of magnitude as the difference bet-ween repeat rune 
-with the same leading edge (see figs. 3 and 6), and can therefore be 
attributed to lack of reproducibility of experimen-bal da-ba. It can be 
concluded that Increasing the leading-edge radius from about 0.0005 to 
about 0.0015 inch has little or no effect cn the measured drag coefficient. 

It is thus reasonable to assume that the discrepancies between 
theoretical and measured drag coefficients found In this inyestlgation 
are not caused by a momentum defect originating in the Immediate 
-vicinity of the leading edge. 

Effect of heat -transfer . - The measured flat-plate drag coefficients 
have been compared -with theoretical drag coefficients for an Insulated 
flat plate. This comparison is not strictly -valid because there -will be 
local heat transfer through the leading edge from the chamfered side of 
the plate to the top of the plate. An estimate based on the theory of 
reference 9 shows, however, that the effect of this heat transfer on skin 
friction is about 0.2 percent and hence may be neglected. The tempera- 
ture along the surface of the plate in tunnel 1 varied by less than 
30 Temperature recovery factors based on this surface temperature had 
a "value of 0.88, which is in good agreement with pre-vioUs eaperiments on 
insulated plates . 
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Effect of prole tip height . - The effect of the total -pres sure prohe 
on measured friction-drag coefficients was investigated hy surveying the 
boundary layer under identical conditions with four probes which had 
outside probe tip heists h ranging from 0.003 to 0.016 inch. It ’tras 
found most convenient to correlate results of these measurements on the 
basis of the difference between measured and theoretical friction-drag 
coefficients . 

This difference 2(0 - %)/z was first plotted against 
which is a function of only the tunnel inlet conditions^ and for given 
values of 2 and h, 2(0 - 6ijt) /x was found to be independent of 

Typical curves of 2(0 - 0g>)/x against are shown in 

figure 18. There is, of course, a large amount of scatter, but it 
should be remembered that an error of about 2 percent in 0 will cause 
an error of up to 30 percent in 2(0 - 0 qi)/x, depending on the magnitude 

of 2(0 - 0^)/x. It is evident, however, that there is no definite 

trend or variation of points in figure 18, and that, in the range of 
Eeynolds number Investigated, 2(0 - 6^)/x can be considered constant 

as long as h and x remain constant. The arithmetic mean of all 
points taken -vriLth a given probe at a given station tTas therefore found, 
and henceforth only this mean value of 2(0 - 0iji)/x will be used. 

A plot of the average friction-drag discrepancy 2(0 - %i)/z 

against probe tip height h is presented in figure 19. A straight line 
is faired through each set of points representing a constant distance x 
from the leading edge. The expression 


2(0 - 0p) 

X 


= E(z) h 


0.64 


(4) 


represents all points reasonably well. Because h is the only variable 
bet^reen points at a given station, eq.uation (4) establishes that the 
greater part, or perhaps all, of the discrepancy between measured and 
theoretical friction-drag coefficients is caused by the finite height of 
the total-pressure probe tip. It is evident, however, that this probe 
effect depends also on the distance from the leading edge. 


In order to determine q.uantltatively this dependency on x, 
measurements were made in tunnel 1 at three values of x with a constant 
probe tip height. A logarithmic plot of the average friction-drag dis- 
crepancy factor against distance from the leading edge is shown in 
figure 20. Points are well represented by a straight line of the form 


2(0 - 


= kx-1-22 


X 


( 5 ) 
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According "bo ©(luatiion (4) ^ "bliis may be generalized bo obber probe beigbbe 
by allowing k bo be a funcbion of b. 

Ib is Inberesbing bo nobe bbab bbe ezponenb of l/z In eq.uabion (5) 
is nearly Urice bbe ezponenb of b in eq.uablan (4) . This facb 
Immediabely suggesbs a correlabion on bbe basis of b/z^. Ib was found 
bbab bbe coefficlenbs of b in eguabion (4) and of z in eq.uabiotn (5) 
are compabible, so bbab for correspcnding values of b and z bbe 
values compubed frcan bobb eq.uabions iirill agree. 

Bobb sebs of daba were bberefore combined and bbe fricbion-drag 
discrepancy facbor 2(0 - ftp)/z 'v/as plobbed cm logaribbmic coordinabes 

againsb b/z^. Tbe sbrai^b line 


2(0 - 0g,) 

z 


0.02 



( 6 ) 


was faired bbrougb all daba polnbs (fig. 21) . 


Tbe polnbs plobbed in figure 21 represenb a range in Reynolds number 

from 10^ bo 10®, a range of z from l/2 bo 3^ inches, and a range in 

8 


probe bip beigbb from 0.003 bo 0.016 incb, and were obbained ab Macb 
numbers of 2 and 3. Eacb poinb in bbis figure corresponds bo bbe 
aribbmeblc mean of from 5 bo 10 measured profiles, so bbab more bban 
60 individual boundary lajnsr surveys are represenbed. Therefore no 
doubb remains bbab bbe fricblon-darag discrepancy facbor as given by 
eq.uabicm (6) is valid for flab plabe fricbion-drag coefficlenbs obbained 
under bbe aforemenbicmed condiblons. 


Ib should be poinbed oub, however, bbab bbe final correlabion was 
based on bbe parameber (b/z2)*^*®^, and bbab bbe coeffioienb of 

(b/z^)®‘®^ in eq.uabion (6) is bberefore nob dimensionless. A correla- 
bion based on bd/z^, where d is bbe t-ridbb of bbe probe opening, would 
be dimensionally correcb and compabible wibb eq.uabion (6) because d 
was nob varied appreciably in bbe presenb besbs. If bbe probe effecb 
is assumed bo be primarily one of Inberference, bben bbe skln-fricbion 
discrepancy facbor would be ezpecbed bo approach zero as d approached 
zero. Bub ib is unlikely bbab a corbelablon based cm bbe probe bip 
area bd could be valid for all bip shapes. Thus, bbe complebe corre- 
labion facbor mlgbb be a funcbicm of bobb bd/z2 and a nondlmensional 
shape facbor. 


A comparison wibb resulbs of references 3 and 12 shows bbab such 
an addiblcmal facbor is req.ulred in bbe correlabion eq.uabion (6) . Tbe 
bhree polnbs from reference 3, plobbed in figure 21, were obbained ab 



16 


KACA TH 2891 


three different values of z with one TOlue of h.. The point from 
reference 12 (also shown in fig. 21) represents tlae arithmetic mean 
of several points obtained at constant h. and x, hut at different 
■values of in-terestlng to note that the friction-drag 

discrepancy factor for the results of reference 12 -was again indepen- 
dent of for Ee^. < 10^. The results of both references are con- 

siderably higher than the friction-drag discrepancy factor found in 
the present study. It appears q.uite possible that a probe tip shape 
factor could bring about agreement for all da’ta, but no such factor 
can be predicted without a systematic ■variation of probe tip width. 

On the other hand^ the fact ■that the leading-edge radius of the 
plate used in reference 3 ■was considerably larger than that of the 
present teste may also acco'unt for the difference in results. 

Of further interest is the fact that if the eaponent in equa- 
tion (6) were 2/3 Instead of 0.64, an eoLuatlon of the following form 
would result; 

2(0 - %) = K (7) 

vdiere K is again dimensional. The factor h/-y/z appears to be 
meaningful because all boundary layer variables increase proportionally 
to ->/z. These same ■variables are also a function of however, 

so that the significance of eq.uation (7) is not clear at psresent. 

A theoretical explanation of the probe effect is beyond the scope 
of the present report. It is reasonable to believe that both a distor- 
tion of the boundary layer caused by the presence of the probe and the 
fact that the probe is in a shear flow rather than a uniform flow 
contribute to the discrepancy. Neither a displacement of the effective 
center of the probe of the magnitude predicted in reference 13 nor an 
area weighting of the pressure field ahead of the probe can account 
for the high measured values of 0. The effect of ■viscosity on probe 
reading (ref. 14) appears to be Important only near the plate surface, 
where the measured poln^ts were replaced by a faired 'curve to the origin, 
so that the momentum thickness remains unaffected. 

The fact that the momentum discrepancy factor is independent of 
U]_/vj^ is in itself difficult to explain. At any given station the 

boundary layer thickness changes as is ■varied. At the same 

time, howe^ver,. the total-pressure gradient in the boundary layer also 
changes. It is not known how these and other factors combine to yield 
a constant ■value of 2(0 - ^)/x. 
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An examination of <i8 and Gj profiles for various prote tip 
heights at constant Eeynolds numher also fails to yield a satisfactory 
explanation of the prohe effect. Typical profiles are shorn in 
figures 22 and 23. The op profiles for prohes 1 and 4 are very 
similar in shape, while the profile for prohe 5 exhibits a shift to the 
ri^t. The peak near the edge of the boundary layer which is present 
in Wo of the profiles becomes more apparent as the probe tip height 
is increased. This peak is also present in the 6j profiles (fig. 23), 

but in general it occurs within the boundary layer. All correlations 
presented herein are based on values of 6j obtained at the value of 
y where ba/bj first vanishes. 

An examination of the factors to which mi^t be attributed the 
discrepancy between measured and theoretical friction-drag. coefficients • 
has therefore showu that only one of these factors, namely, the probe 
effect, leads to a positive result. The fact that profiles measured 
under identical conditions but -vrlth probes of different tip heights 
yield distinctly different values of 6 is evidence enough to establish 
this effect. A correlation as given by eq.uaticn (6) further indicates 
that the entire friction-drag discrepancy may be attributed to the 
probe effect, and that if a hypothetical probe of vanishing tip 
height were employed, theoretical values of 9 would be measured^ 

Further substantiation of the oonclusicn that the fl'iction-drag 
discrepancy may be attributed to the probe effect is obtained from 
reference 5, where the results of direct measurements of local skin- 
friction coefficients using the floating-element techniq,ue are reported. 
These measurements were made at a Mach number of 2.55 on a flat plate 
•VTlth a leeidlng-edge radius of less than 0.0005 inch, local skin- 
friction coefficients, when plotted against a Eeyaolds number based on 
the distance x from the leading edge, agree exactly -vrith the theory 
of Chapman and Eubesln for = 2.55. Average friction-drag coeffi- 
cients obtained by integrating these local values would therefore also 
agree ■v/lth theory. Thus, a friction-drag discrepancy factor of zero is 
obtained when no probe is present in the boundary layer. 


CONCmDING EEMAEiES 

Pressure-probe measurements were made in the laminar boundary layer 
developing on a flat plate at free-stream Mach numbers of 2 and 3. High 
average measured flat-plate friction-drag coefficients, which had 
previously been observed, were again found in the present study. 

Possibilities for error such as the effect of the total-pressure 
probe on boundary layer measurements, heat transfer through the leading- 
edge region, a change in leading-edge radius and the strength of the 
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leading-edge "wave, and the possihllity of early transition to turbulent 
flow or "bursts of turbulence" were Investigated for plates with very 
sharp leading edges. Of these factors, only the probe effect led to 
positive results . 

Total -pressure probes of different tip heights placed in a laminar 
boundary layer developing under identical canditicns were found to 
measure different values of average friction-drag coefficients. An 
empirical coirelatlon of these measurements Indicates that a probe of 
vanishing height would measure the theoretically predicted values of 
average flat plate friction-drag coefficients. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, november 14, 1952 
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APPENDIX A 
SYMBOIS 


The followlag symbols are used in this report: 


Px 


average friction-drag coefficient, — 

2 

local friction-drag coefficient (eq.. (3)) 




pitot pressure 

outside height of probe tip 

constdnts of proportionality 

Mach number 

total pressure 

static pressure 

unperturbed pressure used in error calculations 
gas constant 

local stream Reynolds number based on distance x 


rsidius of curvature of leading edge 
total ten^jerature 
ambient temperature 
velocity in x-direction 

chordwlse distance measured from leading edge parallel to stream 
distance measured perpendicular to plate surface 
dlmensionlees velocity parameter 
dimensionless mass-flow parameter 
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6* 'boundary layer displacement tMckness 

Tj^ temperature recovery factor 

9 "boundary layer momentum thiclmess 

X parameter defined by eq.. (B8) 

V kinematic viscosity 

p mass density 

'z shear stress 

Subscripts: 

1 free-stream conditions 

M measured q^uantity 

s upstream stagnation conditions 

T laminar boundary layer theory 

•w Conditions at plate surface 

y conditions evaluated at distance y from plate surface 
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APPENDIZ B 

! 

EEDBCTION OP DAPA 
By Jack M. lande 


PBe easier imental data were evaluated "by lEM macliliLeB utilizing a 
syatem of data reduction capatle of handling many Q3nsBes of gao flow 
problemB. The application of thiB ByBtem to the reduction of boundary 
layer data 1 b diBcuBBed herein. Althou^ only the caee of the laminar 
boundary layer with constant Btatlc preBBure and constant total temper- 
ature nas conBldered in the main part of the report,. the method of cal- 
culation presented in this appendla: applies to both lam i nar and turbulent 
boundary layers with varying static pressure and total temperature. 

In general, the folloidng q.uantities are measured in boundary layer 
Btudies ; 

H(y) pitot preBsure 

p(y) static pressure 

%(y) measured temperature 

Telocity and temperature profiles, in addition to momentum and displace- 
ment thickness, are desired as a result of calculatians . The calcula- 
tions can be greatly simplified with the aid of the following non.dimen- 
sional quantities (ref. 8) : 


u 

a = — — 

(Bl) 


„ pu-Vm 

(B2) 

op = ^ 
P ■ 

(B3) 


Each of these q.uantities is uniq.uely related to the ratio p/s. 


The momentum and displacement thicknesses for Sp/By = 0 and 
"f 0 become 

I -6- a, . 


e == 




(B4) 
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b* 





(B5) 


Tlie mcmentum thickness 
lation of sk±a friction for 


( 


(op + l)p dy 


is req.uired. 


9 is no longer significant in the calcu- 
Sp/Sy 0; in this case the integral 


The first step in the calculation of the aforementioned quantities 
is the temperature and manometer fluid correction of measured pressures 
and the determination of p/H. 


Calculation of a . - The velocity parameter 
hy the isentropic flovr equation 



a is related to 


p/H 


(B6) 


for Mach numbers less than 1, and hy the Eayleigh aquation 

_ (56a? - 7)^ (7 - 
\Hy " (36)® 


(B7) 


for Mach numbers greater than 1. (The ratio of specific heats y -^ms 
taken as 1.4 in eqs. (B6) and (B7).) 


The difficulties inherent in calculating a from these expressions 
are overcome hy the definition and calculation of an Intermediate 
quantity x which is defined hy 

a = V2(l - p/H) [l + x(l - p/s)] (B8) 

A table deck of x and the slope of x at a given point as a 
function of p/H can he prepared -with the aid of equations (B6) to (B8). 

For a given boundary layer profile, x can then he obtained for each 
value of p/H hy linear Interpolation in the table deck. Finally, a 
is obtained hy substitution of corresponding values of x and p/H into 
equation (B8) . 

Calculation of p . - The quantities a and 6 can he related through 
the equation of state 


p = pEt 


(B9) 
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and the one -dimensional energy eq.uation, 'wliicli can he expressed as 

■fc 7 - 1 2 / ^ 

T = 1 - 27 (BIO) 

For a constant value of specific heat and 7 = 1.4, eq.uations (Bl), 
(B9), and (BIO) yield the following es^ressicn for P; 

P = (BID 

1 - 


This expression was used to calculate p from previously deteimlned 
values of a. 

Calculation of op . - The product of eq.uations (Bl) and (Bll) gives 
for ccp 

op = (B12) 

1 - or 


Calculation of p / - The measured temperature Tf^ is related 
to the local static temperature hy 


7-1 ? 

~ = 1 + ^ (B13) 


where t)^ is the temperature prohe recovery factor. The following 

expression is obtained from eq.uations (B13) and (BIO) and the definitions 
of exp : 


P 

Vt 




(1 + lit 


p). 


% 


(B14) 


Calculation of Integrals . - The integrations of the q.uantities 

P/-\/t, op, and otpp at a distance from the wall are carried out using 
Stirling's parabolic rule for both the laminar and turbulent cases. 

Near the irall, however, where data points are unreliable, two different 
procedures are req.uired for the laminar and turbulent cases. 


Because the velocity profile near the wall for a lamineir boundary 
layer is nearly linear, and because the density under conditions of zero 
heat transfer is nearly constant, the following expressions are used; 
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8 = ay 
o8 = by^ 


> (B15) 


These results, together with the assumption of constant p and T, are 
therefore used to evalixate the Integrals for the first points in a laminar 
boundary layer. 


For turbulent boundary layers the form of 8 and ct8 near the irall 
is determined empirically. The following ea^resslcns, based on available 
turbulent boundary layer data, ai^ suggested for 8: 


or 



+ agy + 



+ a^y^ 


p = yiy^^a^y + agy^ 


+ a- 




&4y 




r (B16) 


The q.uantlty o8 is usually well represented by 


oB = bj^y + bgy^ + b3y^ + b^y^ 


(B17) 


Eq.uations (B16) and (B17) together vrlth the assumptions of constant 
static pressure and constant total temperature are used to evaluate the 
integrals for the points near the wall in the turbulent case. 

Combinatione of these Integrals yield the momentum thickness 9 
and the displacement thickness 8*. 

Significant digits . - All calculations are carried out to five or 
more digits, which is sufficient in relation to the Instrumentation used 
in obtaining the experimental data. 

Calculation time . - The actual machine time req.ulred to obtain 
results- for a given boundary layer profile is minor. It has been found, 
however, that data obtained during a given morning’s irind tunnel run 
(about 4 profiles of about 35 points each) can be reduced by one 
operator during the eiftemoon of the same day. In that manner results 
of a set of tests are knom befoie the next set of tests is begun. Many 
more profiles could be handled at the same time, but for the present 
tests the feature of immediate data analysis outireighed the economy of 
accumulating a great amount of data before reduction -was begun. 
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Figure 3. - Typical velocity parameter profiles. Distance from leading edge, 2.38 inches; 
probe tip height, 0.008 inch; free-stream llaoh number, 3. 
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Figure 6. 


- Typical momentum thickness profiles. Distance from leading edge, 2.38 Inches; 
probe tip height, 0.008 Inch; free-stream Mach number, 3. 
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Figure 8. - Pitot pressure In free stream. Distance from leading 
edge, 2.38 inches; distance frcm plate surface, 0.5 inch; free- 
stream Mach number, 3. 
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Figure 10. - Error in computing mcmentum tlilckaesB caused "by total-temperature 

assumption. 
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Error in prote zero setting, in. 

Figure 11. - Effect of error in proibe zero setting on momentum thick- 
ness. Based on theory of reference 9. Disteince from leading edge, 
2.38 inches; free- stream Mach number, 3; local stream Reynolds 
number, 450,000. 
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Figure 13. - Curves indicating treinsition obtained by varying at 

constant distance frcm leading edge and from plate surface. Distance 
from leading edge, 3.38 Incbesj probe tip height, 0.005 inch; free- 
stream Mach number, 3. 
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Figure 15. - Typical velocity profiles obtained in t unn el 1 at 
0° and -5° angle of attack. Free -stream. Mach number, 2. 
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Figure 17. - Momentum thlcikness profiles vith sharp nni^ dull leading edges. 
Distance from leading edge, 3.38 Inches; free-stream Mach number, 3. 





Figure 16. - Friction-drag dlacrepancy 2(6 - Sgi)/5C a£ function of ui/vi. 
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Figure 21. - Frlctloa-drag dlson 
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I Figure 22. - Kanentum paraiuetar profllas for ttirw dlffarant pro'bes. Dlrbane* fram leading edge, 2.58 incheB^ 

local Btrean HeynolilB nmibor, 0.45X10®) froe-Btream Mach nusiber, 3. 
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